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'The  importance  of  delay  (or,  retardation  in  the  rate  of 
fatigue  crack  growth)  produced  by  load  interactions  in  variable- 
amplitude  loading  on  the  accurate  prediction  of  fatigue  lives  of 
aircraft  and  other  engineering  structures  has  been  recognized  for 
some  time  and  has  begun  to  receive  greater  attention  in  recent . 
years.  Recent  investigations  showed  that  the  effects  of  delay  can 
be  quite  large,  and  that  these  effects  need  to  be  taken  into  account 
in  developing  improved  fatigue  analysis  procedures  for  aircraft  and 
other  engineering  structures. 

A number  of  models  (based  on  the  concepts  of  crack  closure, 
residual  stress  intensity  factor,  etc.)  have  been  proposed  to  ac- 
count for  the  effects  of  delay.  These  models,  while  successfully 
predicting  trends  in  the  rate  of  fatigue  crack  growth  for  ran- 
domized load  spectra,  appear  to  break  down  for  ordered  spectra. 
Several  basic  problems  contributed  to  the  lack  of  complete  success, 
and  needs  to  be  resolved  in  the  development  of  improved  models  for 
predicting  load  interaction  effects  (chiefly  delay)  on  fatigue 
crack  growth.  They  are; 
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o Proper  characterization  and  physical  understanding 
of  the  complex  phenomena  of  crack  acceleration  and 
delay  associated  with  changes  in  load  level,  and  of 
the  effects  of  thermal  and  chemical  environments  on 
these  phenomena. 

o Adequate  description  of  crack  tip  stress  intensity 
factors  to  account  for  various  types  of  loading, 
crack  geometry  and  residual  stresses. 

o Proper  description  of  the  kinetics  of  fatigue  crack  growth 
under  constant-amplitude  loading,  including  the  effects 
of  stress  ratio,  cyclic  load  frequency  and  service 
environment,  and  assessment  of  variability  in  these 
kinetic  data. 

To  address  some  of  the  issues,  the  following  three  inves- 
tigations were  carried  out  under  this  grant; 

(1)  Experimental  evaluation  of  crack  closure  and  its 
viability  as  a model  for  describing  fatigue  crack 
growth  kinetics  and  delay. 

(2)  Excuninations  of  the  influences  of  plate  thickness, 

K level  and  chemical  environment  on  fatigue  crack 
growth  response  following  a single  high-load  excur- 
sion (overload)  to  develop  further  phenomenological 
understanding  of  load  interaction  effects  in  fatigue. 

(3)  Assessment  of  the  contributions  of  crack  length 
measurement  interval  and  precision,  and  data  process- 
ing procedures  on  variability  to  provide  background 
for  interpretation  of  experimental  data  on  fatigue 
crack  growth  kinetics  and  their  utilization  in  life 
prediction. 

2219-T851  aluminimi  alloy  plates  were  used  in  the  first  two 
investigations.  Computer  simulation  was  used  in  the  last  inves- 
tigation to  preclude  all  other  variables.  Results  from  the  first 
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investigation  have  been  included  in  an  Interim  Report,  entitled 
"Crack  Closure  and  Fatigue  Crack  Growth  in  2219-T851  Aliaminum 
Alloy",  by  K.  D.  Unangst,  T.  T.  Shih  and  R.  P.  Wei,  dated  August 
1976  (AFOSR-TR-76-1247,  ADA633082) . The  principal  results  have 
been  reported  also  in  a technical  publication,  entitled  "Crack 
Closure  in  2219-T851  Aluminum  Alloy",  by  the  same  authors  in  the 
Journal  of  Engineering  Fracture  Mechanics  (v.  9,  1977,  pp.  725-734). 
Results  from  the  other  two  investigations  have  been  prepared  in 
the  form  of  technical  publications  and  are  included  here  as  Appen- 
dices A and  B.  Brief  summaries  of  these  investigations  are 
given  here. 

The  effects  of  specimen  thickness,  stress  ratio  (R)  and 

maximum  stress  intensity  factor  (K  on  crack  closure  (or  open- 

max 

ing)  and  on  fatigue  crack  growth  kinetics  were  studied  using  a 
2219-T851  aluminum  alloy.  The  crack  length  and  the  occurrence 
of  crack  closure  were  measured  by  an  electrical  potential  method. 

The  experimental  work  was  carried  out  within  the  fraunework  of 
linear-elastic  fracture  mechanics.  The  experimental  results  show 
that  the  onset  of  crack  closure  (or  opening)  depends  on  R, 
and  specimen  thickness.  In  terms  of  the  "effective  stress  in- 
tensity range  ratio"  (U) , as  defined  by  Elber,  the  results  show 
that  U tends  to  increase  for  increasing  R,  decrease  for  increasing 
^max'  decrease  with  increasing  specimen  thickness.  Prom  these 
trends,  it  is  shown  that  the  "effective  stress  intensity  range" 
(AKeff)  does  not  always  increase  with  increasing  stress  intensity 
range  (AK) . The  fatigue  crack  growth  data  show  that  the  specimen 
thickness  does  not  have  a significant  effect  on  crack  growth  in 
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this  material  over  the  lower  crack  growth  rate  region;  below 
about  5 X lO”®  in. /cycle  (1.3  x 10~^)  cm/cycle) . In  the  higher 
crack  growth  rate  region,  above  about  5 x lO”^  in. /cycle  (1.3  x 
lO”^  cm/cycle) , the  crack  growth  rates  are  higher  for  the  thick- 
er specimens.  The  viability  of  the  crack  closure  model  is  ques- 
tioned. The  experimenal  results  show  that  crack  closure  cannot 
, fully  account  for  the  effects  of  stress  ratio  and  specimen  thick- 
ness or  K on  fatigue  crack  growth.  The  use  of  AK  ..  as  a para- 
max  eft 

meter  for  characterizing  the  mechanical  driving  force  for  fatigue 
crack  growth  is  questioned. 

To  develop  further  phenomenological  understanding  of  load 
interaction  effects  in  fatigue,  examinations  of  the  influences 
of  plate  thickness,  stress  intensity  (K)  level  and  chemical  en- 
vironment on  fatigue  crack  growth  response  following  a single 
high-load  excursion  (overload)  were  carried  out  on  a 2219-T851 
aluminum  alloy.  An  overload  ratio  (that  is,  the  ratio  between 
the  magnitude  of  the  overload  and  the  maximum  in  the  subsequent 
constant-cunplitude  fatigue  loading)  of  2.0  was  used.  Experiments 
were  carried  out  in  dehumidified  argon,  air,  and  3.5  pet  NaCl 
solution  at  room  temperature.  The  results  showed  that  delay  (as 
measured  by  the  duration  of  overload  affected  crack  growth)  in- 
creased with  decreasing  plate  thickness  and  with  increasing  K 
level,  and  decreased  with  increasing  aggressiveness  of  the  chemical 
environment. 
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The  high-load  excursion  (overload)  affected  crack  growth 
through  a region  of  material  ahead  of  the  crack  tip.  Within  this 
overload  affected  zone,  crack  growth  rate  first  increased  (some- 
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times) , followed  by  fairly  rapid  decrease  to  a minimum  value 
(delayed  retardation) , and  then  increased  gradually  to  its  steady- 
state  value.  The  overload  affected  zone  size  was  found  to  depend 
on  crack-tip  constraint  and  to  be  independent  of  chemical  envi- 
ronment, and  was  found  to  be  equal  to  the  appropriate  (plane- 
strain  or  plane-stress)  plastic  zone  size  for  the  overload.  Iden- 
tification of  a delayed  retardation  zone  was  made,  and  identifi- 
cation of  this  zone  with  the  cyclic  plastic  zone  size  for  the  pre- 
ceding fatigue  loading  was  suggested.  The  effects  of  plate  thick- 
ness, K level  and  chemical  environment  on  delay  were  considered  in 
relation  to  their  respective  influences  on  the  overload-affected- 
zone  and  delayed-re tarda tion- zone  sizes,  and  on  the  rate  of  fa- 
tigue crack  growth.  A residual  stress  intensity  concept  for  des- 
cribing fatigue  crack  growth  response  within  the  overload  affect- 
ed zone  was  considered.  With  suitable  modifications,  reasonable 
estimates  of  crack  growth  response  could  be  obtained.  Further 
verification  and  understanding  of  these  modifications  are  discuss- 
ed. 

A computer  simulation  of  the  primary  fatigue  crack  growth 
data  (crack  length  versus  elapsed  cycles  or  a versus  N)  was  made 
to  examine  the  influences  of  crack  length  measurement  interval 
and  precision,  and  of  the  associated  data  processing  procedure  on 
the  resulting  data  on  fatigue  crack  growth  kinetics  (i.e.,  da/dN 
versus  AK) . Variability  in  the  derived  growth  rate  (da/dN)  data 
depended  strongly  on  the  magnitude  of  the  measurement  interval 
relative  to  the  measurement  precision.  It  was  reduced  by  those 
data  processing  procedures  (such  as  the  incremental  polynomial 


methods)  that  fitted  a smooth  curve  through  portions  of  the  pri- 
mary data.  Such  procedures,  however,  introduced  significant  bias 
into  the  derived  data  on  fatigue  crack  growth  kinetics.  The  re- 
sults further  suggested  that  much  of  the  variability  in  the  pub- 
lished data  on  fatigue  crack  growth  kinetics  might  be  attributed 
to  sources  described  herein.  Caution  is,  therefore,  recommended 
in  attempting  to  draw  statistical  inferences  regarding  material 
variability  from  these  data. 

Results  from  these  investigations  have  provided  additional 
insignt  into  load  and  environment  interactions  in  fatigue  crack 
growth  under  variable-cimplitude  (spectrum)  loading,  as  well  as 
for  fatigue  under  constant-amplitude  loading.  These  results 
should  be  considered  in  assessing  current  and  improved  life  pre- 
diction procedures. 
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ABSTRACT 

To  develop  further  phenomenological  understanding  of  load 
interaction  effects  in  fatigue,  examinations  of  the  influences 
of  plate  thickness,  stress  intensity  (K)  level  and  chemical  en- 
vironment on  fatigue  crack  growth  response  following  a single 
high-load  excursion  (overload)  were  carried  out  on  a 2219-T851 
aluminum  alloy.  An  overload  ratio  (that  is,  the  ratio  between 
the  magnitude  of  the  overload  and  the  maximum  in  the  subsequent 
constant-amplitude  fatigue  loading)  of  2.0  was  used.  Experiments 
were  carried  out  in  dehumidified  argon,  air,  and  3.5  pet  NaCl  so- 
lution at  room  temperature.  The  results  showed  that  delay  (as 
measured  by  the  duration  of  overload  affected  crack  growth)  in- 
creased with  decreasing  plate  thickness  and  with  increasing  K level, 
and  decreased  with  increasing  aggressiveness  of  the  chemical  envi- 
ronment. 
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The  high-load  excursion  (overload)  affected  crack  growth 
through  a region  of  material  ahead  of  the  crack  tip.  Within  this 
overload  affected  zone,  crack  growth  rate  first  increased  (some- 
times) , followed  by  fairly  rapid  decrease  to  a minimum  value 
(delayed  retardation) , and  then  increased  gradually  to  its  steady- 
state  value.  The  overload  affected  zone  size  was  found  to  de- 
pend on  crack-tip  constraint  and  to  be  independent  of  chemical 
environment,  and  was  found  to  be  equal  to  the  appropriate  (plane- 
strain  or  plane-stress)  plastic  zone  size  for  the  overload.  Iden- 
tification of  a delayed  retardation  zone  was  made,  and  identifi- 
cation of  this  zone  with  the  cyclic  plastic  zone  size  for  the 
preceding  fatigue  loading  was  suggested.  Tl  • effects  of  plate 
thickness,  K level  and  chemical  environment  n delay  were  con- 
sidered in  relation  to  their  respective  influences  on  the  over- 
load-affected-zone and  delayed-retardation-zone  sizes,  and  on  the 
rate  of  fatigue  crack  growth.  A residual  stress  intensity  concept 
for  describing  fatigue  crack  growth  response  within  the  overload 
affected  zone  was  considered.  With  suitable  modifications,  rea- 
sonable estimates  of  crack  growth  response  could  be  obtained. 
Further  verification  and  understanding  of  these  modifications  are 


discussed. 


INTRODUCTION 


The  importance  of  delay  (retardation  in  the  rate  of  fatigue 
crack  growth) , produced  by  load  interaction  in  variable-2unplitude 
loading,  on  the  accurate  prediction  of  fatigue  lives  of  aircraft 
and  other  engineering  structures  has  been  well  recognized  for  some 
time  [1-6] . Several  models  have  been  proposed  to  account  for  the 
effects  of  delay  in  fatigue  analysis  of  structures  [7-15].  These 
models,  while  successfully  predicting  trends  in  the  rate  of  fatigue 
crack  growth  for  randomized  load  spectra,  appear  to  break  down 
(in  most  cases)  for  ordered  spectra  [6-18] . Development  of  im- 
proved models  to  account  for  the  effects  of  delay,  however,  depends 
on  a more  complete  understanding  of  the  process  of  delay  itself. 
Available  systematic  data  for  simple  loading  conditions  (particu- 
larly for  the  Ti-6A1-4V  alloy)  indicate  that  delay  in  fatigue 
crack  growth  is  a highly  complex  phenomenon,  and  can  be  affected 
by  both  chemical  and  thermal  environments,  as  well  as  by  section 
size  (specifically  plate  thickness)  [19-26]  . To  expand  on  the  phe- 
monenological  understanding  of  delay,  the  effects  of  chemical  en- 
vironment, (dehumidified  argon,  air  and  3.5  pet  NaCl  solution) 
'tress  intensity  (K)  level  and  plate  thickness  were  investigated 
for  fatigue  crack  growth  following  a single  high-load  excursion 
("overload")  in  a 2219-T851  aluminum  alloy  at  room  temperature. 

Delay  properly  refers  to  the  period  of  abnormally  low  rate, 
or  approximately  zero  rate,  of  fatigue  crack  growth  between  a de- 
crease in  load  level  and  the  establishment  of  a rate  of  crack  grow- 
th commensurate  with  that  for  constant-aunplitude  loading  at  the 
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prevailing  (lower)  load;  that  is,  between  Points  B and  E in  the 
schematic  diagram  shown  in  Figure  1.  (Note  that  in  Figure  1 
the  phenomena  of  crack  acceleration  and  delayed  retardation,  that 
is,  B to  B'  and  B'  to  C,  are  also  illustrated  [19]).  The  response 
following  a single  high-load  excursion  (overload)  is  depicted 
here  (see  inset  in  Figure  1 where  the  loading  parameter  is  given 
in  terms  of  the  crack-tip  stress  intensity  factor  K) . Following 
Jon^s  and  Wei  [21] , subscripts  1 and  2 are  used  to  denote  parame- 
ters associated  with  the  overload  and  the  steady-state  fatigue 
loads  respectively. 

Two  of  the  current  procedures  used  to  define  delay  in  terms 
of  the  number  of  elapsed  cycles  are  illustrated  in  Figure  la.  In 
the  first  method,  delay  (N^)  is  artificially  defined  as  a period 
of  zero  crack  growth,  represented  by  the  interval  B'E'  in  Figure 
la,  and  is  obtained  by  extrapolating  the  constant-amplitude  grow- 
th curve  (Curve  EF)  to  E' . This  definition  of  delay  was  used  with 
considerable  success  for  the  Ti-6A1-4V  alloy,  since  there  was  min- 
imal crack  acceleration  and  delayed  retardation  [19] . The  alter- 
native definition  for  delay  includes  the  entire  period  affected 
by  the  high-load  excursion,  and  delay  (N*)  is  depicted  by  the 
interval  from  B to  E (see  Figure  la) . A corresponding  "overload 
affected  zone"  is  defined  in  Figure  lb  (Point  B'  to  E)  along  with 
the  detailed  crack  growth  response. 

Because  of  the  occurrence  of  crack  acceleration  and  delayed 


retardation  in  some  of  the  aluminum  alloys  and  to  provide  more 
detailed  description  of  crack  growth  response,  the  second  of 
these  two  procedures  was  implicitly  adopted  for  this  investigation. 


Emphasis  was  placed,  therefore,  on  examining  the  rate  of  fatigue 
crack  growth  as  a function  of  crack  extension  and  of  elapsed 
cycles  before  and  following  a single  high-load  excursion  (over- 
load) , see  Figure  1.  Experiments  were  carried  out  within  the 
fraunework  of  linear  fracture  mechanics.  The  crack  tip  stress  in- 
tensity factor  K or  AK  was  used  to  characterize  the  mechanical 
driving  force  for  crack  growth. 

MATERIAL  AND  EXPERIMENTAL  WORK 

Material 

A 7.6  cm  (3  in.)  thick  plate and  a 1.65  cm  (0.65  in.) 

(2) 

thick  plate  of  2219-T851  aluminum  alloy  were  used  in  this 

study.  Chemical  composition  and  tensile  properties  for  the  7.6 

(3) 

cm  thick  plate  are  given  in  Table  1 . 

Specimen  and  K Calibration 

Wedge-opening- load  (WOL)  specimens,  with  half-height  to  wid- 
th ratio  (H/W)  of  0.486,  were  selected  for  use  in  this  study. 

The  specimens  were  oriented  in  the  longitudinal  (LT)  orientation. 
Two  slightly  different  configurations  were  used  to  facilitate 
attachment  of  leads  for  the  electrical-potential  crack  measurement 
system  (see  Figure  2).  Specimens  with  thicknesses  of  0.25,  0.51, 
0.76  and  1.27  cm  (0.1,  0.2,  0.3  and  0.5  in.)  were  machined  from 


(1)  Plate  supplied  by  Westinghouse  Electric  Corporation,  and  was 
a part  of  the  material  used  in  an  AFML  sponsored  program 
(Contract  F33615-75-C-5064)  on  fatigue  crack  growth. 

(2)  Plate  furnished  by  AFFDL. 

(3)  Chemical  composition  and  tensile  properties  were  not  deter- 
mined on  the  1.65  cm  (0.65  in.)  thick  plate.  They  are  ex- 
pected, however,  to  be  comparable  to  those  for  the  7.6  cm 
(3.0  in.)  thick  plate. 
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•the  7.6  cm  (3  in.)  thick  plate.  These  specimens  conformed  to 
dimensions  shown  in  Figure  2a,  and  were  randomized  in  the  thick- 
ness direction.  Specimens  from  the  1.65  cm  (0.65  in.)  thick 
plate  were  machined  in  accordance  with  Figure  2b,  and  were  in 
full  thickness.  An  initial  (or  crack  starter)  notch,  1.96  cm 
(0.77  in.)  in  length,  was  introduced  into  each  specimen  by  elec- 
tro-discharge machining  (EDM) . Each  specimen  was  precracked  in 
fatigue  through  a decreasing  sequence  of  loads  that  terminated 
at  the  desired  load  level  (or  initial  K)  for  the  actual  experiment. 
The  precracking  procedure  provided  a fatigue  crack  of  about  0.33 
cm  (0.13  in.)  in  length  from  the  starter  notch,  corresponding 
to  a crack  length  of  about  2.29  cm  (0.9  in.)  at  the  start  of  each 
experiment.  This  precracking  procedure  ensured  that  the  subse- 
quent fatigue  crack  growth  would  be  through  material  that  had  not 
been  altered  by  the  notch  preparation  procedure,  and  would  be  un- 
affected by  the  starter  notch  geometry. 

Stress  intensity  factor,  K,  for  the  WOL  specimen  was  computed 
from  Eq.  1 [27,28]: 


K = ^/a  (30.96  - 195.8  (a/W)  + 730.6  (a/W)  ^ 
- 1186.3  (a/W)^  + 754.6  (a/W)^] 


(1) 


P * applied  load;  B = specimen  thickness;  W = specimen  width;  and 
a * crack  length.  Both  specimen  width  and  crack  length  were  meas- 
ured from  the  line  of  loading,  as  shown  in  Figure  2. 


Experimental  Procedures 

Fatigue  crack  growth  experiments  to  study  load  interaction 


i 
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effects  were  carried  out  in  a closed-loop  electrohydraulic  test- 
ing machine  operated  in  load  control.  Load  control  was  estimated 
to  be  better  than  ± 1 pet.  The  effects  of  plate  thickness,  stress 
intensity  (K)  level  and  chemical  environment  on  delay  (or  fatigue 
crack  growth  response) , at  room  temperature,  following  a single 
high-load  excursion  or  overload  (see  inset  in  Figure  1)  were  exa- 
mined. The  investigation  was  limited  to  an  overload  ratio 
^2max^  of  2.0,  with  the  intervening  fatigue  loading  at  a load 
ratio  (R)  of  0.05. 

Two  different  procedures  were  used  to  accommodate  a number 
of  overload  tests  on  each  specimen.  In  the  first  procedure,  an 

attempt  was  made  to  achieve  the  seune  K,  and  K.  at  each  over- 

^ Imax  2max 

load  point,  and  shall  be  designated  hereafter  as  the  "constant  K" 
test.  In  the  second  procedure,  the  magnitudes  of  the  loads  were 
maintained  constant;  that  is,  "constant- load"  test.  In  the  con- 
stant K tests,  the  fatigue  load  magnitudes  were  selected  to  attain 
a desired  value  of  at  a pre-selected  crack  length.  Fatigue 

loading  was  then  interrupted,  and  a single  overload  cycle  was 
applied,  using  the  single  cycle  feature  of  the  testing  machine, 
with  a triangular  waveform  at  0.01  Hz.  Following  the  overload, 
fatigue  loading  was  resumed  at  the  Scime  loading  conditions  immedi- 
ately before  the  overload  application  and  was  continued  to  extend 
the  crack  beyond  the  point  where  steady-state  rates  were  re-es- 
tablished (that  is,  beyond  Point  E in  Figure  1).  At  this  point, 
the  fatigue  loads  were  reduced  such  that  the  desired  would 

be  achieved  at  the  next  overload  point.  The  change  in  K level 
within  any  overload-affected-zone  was  estimated  to  be  less  than 


10  pet.  The  procedure  was  repeated  until  the  desired  number  of 
overload  tests  were  accomplished,  and  provided  a number  of  tests 


at  the  Scune  K2max  constant  load  test,  the  fatigue 

load  magnitude  was  kept  the  scune  for  each  specimen.  Overloads 
were  applied  in  the  saune  manner  as  that  for  the  constant  K tests. 
With  this  procedure,  a number  of  values  were  examined  using 

a single  specimen. 

Both  of  the  procedures  were  used,  in  conjunction  with  a dc 
crack  monitoring  system,  in  studies  on  the  7.6  cm  (3.0  in.)  thick 
plate,  and  involved  manual  interruption  and  reinitiation  for  each 
of  the  overload  tests.  Fatigue  cracking  between  overload  was 
carried  out  under  sinusoidally  varying  loads  at  10  Hz.  For  the 
1.65  cm  (0.65  in.)  thick  plate,  only  the  constant  K procedure  was 
used.  Crack  growth  was  measured  by  an  ac  crack  monitoring  system. 
The  tests  were  fully  automated  and  were  controlled  by  an  inter- 
active digital  computer.  Fatigue  cracking  between  overloads  was 
carried  out  under  sinusoidally  varying  loads  at  20  Hz,  and  the 
overloads  were  applied  using  a triangular  waveform  at  0.01  Hz. 


Crack  Monitoring  System 

Both  a dc  and  an  ac  electrical  potential  system  were  used 
for  monitoring  crack  growth  [29-32] . The  dc  system  was  used  with 
specimens  from  the  7.6  cm  (3.0  in.)  thick  plate,  and  the  ac  sys- 
tem was  used  for  specimens  from  the  1.65  cm  (0.65  in.)  thick  plate. 

For  the  specimen  geometry,  analytical  relationships  between 
crack  length  and  electrical  potential  were  not  available,  and  ex- 
perimental calibration  curves  had  to  be  established.  Experimental 
calibrations  were  accomplished  either  by  making  simultaneous 
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visual  and  electrical  potential  measurements  of  crack  length  on 
specimens  fatigued  in  air,  or  by  comparing  electrical  potential 
measurements  against  fatigue  markings  on  fracture  surfaces  pro- 
duced by  intentionally  introduced  high-load  excursions  (overloads) 
during  fatigue  in  air.  Calibration  results  are  given  in  Figures 
3 and  4. 

The  results  in  Figure  3a  are  dc  measurements  for  the  specimen 
configuration  shown  in  Figure  2a  at  three  specimen  thicknesses, 
and  are  given  as  crack  length  (a)  versus  normalized  potential 
(V*)t  These  results  show  the  reproducibility  between  specimens 
and  confirm  that  the  calibration  curve  is  independent  of  specimen 
thickness.  The  following  second  degree  polynomial,  Eq.  2,  pro- 
vided the  best  (least-square)  fit  to  the  data,  and  was  used  as 
the  calibration  curve; 

a = 0.792  + 3.43V*  - 1.54V*^  (a  in  in.) 

5 (2) 

a = 2.01  + 8.71V*  - 3.91V*'‘  (a  in  cm) 

where  a = crack  length,  V*  = (V-V^)/V^,  V^  = reference  potential 
associated  with  the  initial  notch,  and  V = V(a)  = potential  at  a 


tThe  electrical  potential  method  provides  measurements  of  crack 
length  averaged  through  the  thickness,  while  the  visual  method 
gives  measurements  of  the  crack  length  at  the  specimen  surface 
only.  Crack  length  measurements  made  by  these  two  methods  would 
differ  because  of  crack  front  curvature.  The  discrepancy  was 
significant  for  the  thicker  specimen.  Corrections  for  crack 
front  curvature  were  made  by  measuring  average  crack  lengths  from 
the  fatigue  markings  (introduced  during  the  calibration  tests  by 
changing  the  load  amplitude)  after  specimen  fracture.  The  aver- 
age crack  length  was  computed  on  the  basis  of  five  measurements  - 
one  at  each  specimen  surface,  one  along  each  of  the  quarter- 
thickness  planes  and  one  along  the  mid-thickness  plane.  The 
"corrected"  crack  lengths  are  used  in  Figures  3 and  4 and  in 
deriving  Eqs.  2,  3 and  4. 
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crack  length  a.  The  calibration  result  in  Figure  3b  is  for  the 
specimen  configuration  shown  in  Figure  2b,  at  1.65  cm  (0.65  in.) 
thickness,  using  the  dc  measurement  system.  The  corresponding 
calibration  curve  is  given  by  Eq.  (3) . 

a = 0.77  + 2.61V*^  - 0.26V*^  - 0.81V*^  (a  in  in.) 

2 3 

a = 1.96  + 6.63V*  - 0.66V*  - 2.06V*^  (a  in  cm) 


Accuracy  of  crack  length  measurement  was  estimated  to  be  better 
than  1 pet  for  crack  lengths  from  about  2 to  4.8  cm  (0.8  to  1.85 
in.).  The  resolution,  based  on  a fixed  working  dc  current  of 
about  10  amperes,  however,  depended  on  specimen  thickness,  and 
was  only  slightly  dependent  on  crack  length.  For  the  1.27  cm 
(0.5  in.)  thick  specimen,  crack  length  resolution  was  better  than 
0.005  cm  (0.002  in.)  based  on  O.lyV  resolution  in  electrical  po- 
tential. Resolution  for  the  other  specimens  varied  in  inverse 
proportion  to  the  specimen  thickness,  that  is,  0.007  cm,  0.002  cm 
and  0.001  cm  (0.003  in.,  0.0008  in.  and  0.0004  in.)  for  the  1.65, 
0.51,  and  0.25  cm  (0.65,  0.2,  and  0.1  in.)  thick  specimens  respec- 
tively. 

Calibration  results  obtained  from  ac  measurements  on  1.65  cm 
(0.65  in.)  thick  specimens,  with  geometry  shown  in  Figure  2b,  are 
given  in  Figure  4 as  crack  length  (a)  versus  potential  difference 
(V-V^) . The  calibration  results  can  be  represented  by  Eq.  (4). 
a = 0.77  + 0.1586(V  - V^)  (a  in  in.) 


(4) 

a * 1.96  + 0.4028(V  - V^)  (a  in  cm) 

V and  V^  are  defined  in  the  same  manner  as  that  of  the  dc  system 
and  are  given  in  microvolts.  Because  of  variations  in  the 
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inductive  contribution  to  V and  from  specimen  to  specimen / 
normalization  of  electrical  potential  could  not  be  used  in  this 
case.  Unlike  Eqs.  (2)  and  (3),  therefore,  Eq.  (4)  is  restricted 
to  the  particular  specimen  and  operating  current.  Accuracy  of 
crack  length  measurement  with  the  ac  system  was  estimated  to  be 
better  than  1 pet  also,  for  crack  lengths  from  about  2 to  4.8  cm 
(0.8  to  1.85  in.).  The  resolution  was  better  than  0.005  cm 
(0.002  in.)  based  on  12.5  nV  resolution  in  electrical  potential. 
Because  of  the  use  of  a "lock-in"  amplifier,  better  signal-to- 
noise  ratio  was  obtained  with  the  ac  system. 

Test  Environment 

The  load  interaction  experiments  were  carried  out  in  dehu- 
midified argon,  air  and  3.5  pet  sodium  chloride  solution  (3.5  pet 
NaCl)  at  room  temperature.  For  tests  in  dehumidified  argon,  the 
environment  was  maintained  around  the  crack  by  flowing  argon 
(purified  by  a suitable  purification  system)  through  cheunbers 
clamped  to  the  faces  of  the  specimen  [31] . Tests  in  air  were  car- 
ried out  with  no  special  control  (with  relative  humidity  in  the 
range  of  30  to  70  pet).  3.5  pet  NaCl  solution  was  obtained  by 
dissolving  NaCl  in  triply  distilled  water.  During  testing,  the 
solution  was  maintained  around  the  crack  in  (inert)  chambers  at- 
tached to  the  specimen. 

EXPERIMENTAL  RESULTS 

Fatigue  crack  growth  response  following  a single  high-load 
excursion  (overload)  is  illustrated  schematically  in  Figure  1. 

The  actual  experimental  data  (electrical  potential  versus  time) 


generally  conformed  to  that  depicted  in  Figure  la,  although  the 
records  invariably  exhibited  perturbations  during  and  immediately 
following  overload  applications*.  Fatigue  crack  growth  response 
data  are  reported  as  the  rate  of  crack  growth  (da/dN)  versus  rela- 
tive crack  position  (Aa..  = a - a.^ ) with  respect  to  the  point  of 
overload  application  (that  is,  in  a form  similar  to  Figure  lb) , 
and  correspondingly  as  da/dN  versus  ANqj^  ^ 

subscripts  "OL"  denote  crack  length  and  number  of  cycles  at  which 
an  overload  was  applied. 

To  facilitate  discussion,  a number  of  formally  defined  plastic 
zone  dimensions  were  computed  and  indicated  on  the  appropriate 
figures.  The  dimensions  include  [33-35] : 

Monotonic  Overload  Plastic  Zone  Size: 


1 Imax 


■pi  IT  a. 


_ _1_  Imax 

•Ipl  3Tr  Cyg 


for  plane-stress 


for  plane-strain 


Monotonic  Plastic  Zone  Size  Associated  With  Fatigue: 


1 2max 


■p2  TT  a 


_L  2max 
■IP2  Bit  a„„ 


for  plane-stress 


for  plane-strain 


*These  perturbations  introduced  uncertainties  into  crack  growth 
and  growth  rate  measurements  immediately  following  the  overloads. 
These  uncertainties  are  discussed  in  the  Discussion  section. 
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Cyclic  Plastic  Zone  Size: 

1 I ^^2  I ^ 

r. 


pc 


1 

IT 


2o. 


ysj 


f AK2 1 ^ 


Ipc  3it 


2 a 


ysJ 


for  plane-stress  (9) 


for  plane-strain  (10) 


a in  Eqs.  (5)  to  (10)  is  the  yield  strength  in  tension,  and 
ys 

the  zone  sizes  are  equal  to  twice  the  formally  calculated  plastic 
zone  correction  factors  (for  ex5unple,  r^j^  * 2 [33-35]  . 

For  convenience,  the  results  for  specimens  from  the  7.6  cm 
(3.0  in.)  thicJc  and  the  1.65  cm  (0.65  in.)  thicJc  plates  are  re- 
ported separately.  Note  that  experimental  wor)c  on  the  7.6  cm 
thiclc  plate  involved  the  use  of  a dc  crac)c  monitoring  system,  and 
required  manual  interruptions  of  the  testing  machine  for  overload 
applications.  Tests  on  the  1.65  cm  (0.65  in.)  thiclc  plate  on  the 
other  hand  were  fully  automated,  and  utilized  an  ac  crack  measure- 
ment system. 


7.6  cm  (3.0  in.)  Thick  Plate 

Fatigue  crack  growth  responses  associated  with  a single 
high- load  excursion  (overload)  , with  an  overload  ratio 
^2max^  of  2.0,  for  specimens  of  different  thicknesses  machined 
from  the  7.6  cm  thick  plate  and  tested  in  different  environments 
are  shown  in  Figures  5 through  10.  A number  of  overload  experi- 
ments were  carried  out  on  each  specimen.  Individual  runs  on  each 
specimen  are  designated  sequentially  by  letters  (A  to  E)  or  by 
numbers  (1  to  5),  and  corresponded  to  increasing  crack  lengths  at 
which  the  overloads  were  applied  • Fatigue  cracking  between 
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overloads  were  carried  out  under  sinusoidally  varying  loads  at 
a load  ratio  of  0.05,  and  frequency  of  10  Hz. 

Data  from  0.25  cm  (0.1  in.)  and  1.27  cm  (0.5  in.)  thicJc 

specimens,  tested  in  dehumidified  argon  using  the  "constant  K" 

-3/2  1/2 

procedure  at  K2max  ~ MN-m  ' (10  ksi-in.  ),  are  shown  in 

Figures  5 and  6.  Two  specimens  were  tested  at  each  specimen  thick- 
ness. There  was  considerable  scatter  in  the  data,  but  the  overall 
trends  were  consistent,  and  did  not  suggest  any  significant  in- 
fluence of  a^j^  (that  is,  from  proximity  of  the  crack  tip  to  the 
back  edge  of  the  specimen) . Two  of  the  specimens  (F7-33-1  and  CO- 
24-5)  exhibited  significant  crack  acceleration,  while  the  other 
two  did  not.  The  growth  rates  decreased  following  the  overload 
(and  crack  acceleration) , and  appeared  to  reach  a minimum  value 
after  growing  a distance  equal  to  about  ^^^2  ^p2  the  over- 

load point.  Thereafter,  the  rates  increased  until  they  re-estab- 
lished their  equilibrium  value.  The  overload  affected  zone  appear- 
ed to  be  equal  to  about  1 rp^  for  the  thinner  specimens,  and  to 
about  2 th®  thicker  specimen  (Figure  5)  . Overall  delay 

appeared  to  be  longer  for  the  thinner  specimen  than  the  thicker 
specimen  (Figure  6)  and  is  consistent  with  previous  results  on  a 
7075-T6  aluminum  alloy  [25].  These  and  other  aspects  of  delay  are 
discussed  in  a later  section. 

Comparable  data  from  0.51  cm  (0.2  in.)  and  1.27  cm  (0.5  in.) 
thick  specimens,  tested  in  dehumidified  argon  using  the  "constant 
load"  procedure,  are  shown  in  Figures  7 and  8.  These  results  are 
consistent  with  those  shown  in  Figures  5 and  6 with  respect  to  the 
effects  of  specimen  thickness  on  the  sizes  of  delayed  retardation 
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region  and  overload  affected  zone,  and  on  delay.  Both  of  the 
zone  sizes  increased  with  increasing 

Overload  affected  fatigue  crack  growth  response  data  for  a 
0.51  cm  (0.2  in.)  thick  specimen,  tested  in  air  by  the  "constant 
load"  procedure,  are  shown  in  Figure  9.  Comparison  of  Figures 
7a  and  9a  shows  that  the  rate  of  crack  growth  (both  inside  and 
away  from  the  overload  affected  region)  was  faster  in  air  than  in 
dehumidified  argon.  The  size  of  the  overload  affected  zone,  how- 
ever, does  not  appear  to  be  affected  by  the  environment  and  is 
equal  to  about  1 The  size  of  the  delay  retardation  region, 

on  the  other  hand,  appears  to  be  much  smaller  for  the  specimen 
tested  in  air  (being  approximately  equal  to  the  cyclic  plastic 
zone  size  t and  delay  appears  to  be  shorter  in  air  (compare 

Figures  7a,  8a  and  9) . The  dependence  of  these  zone  sizes  on  K 
level  are  reflected  through  their  dependence  on  the  plastic  zone 
sizes  (rpj^  and  r^p^)  . 

The  effects  of  environment  on  fatigue  crack  growth  response 

in  this  material  are  more  readily  seen  in  Figure  10.  Data  from 

0.51  cm  (0.2  in.)  thick  specimens  tested  in  dehumidified  argon  and 

in  air,  and  for  0.76  cm  (0.3  in.)  thick  specimen  tested  in  3.5  pet 

-3/2  1/2 

NaCl  solution  are  compared  at  = 10.8  MN-m  ' (9.8  ksi-ii.  ' ). 

The  results  again  suggest  that  the  overload  affected  zone  was 
essentially  independent  of  environment.  The  rates  of  crack  grow- 
th, on  the  other  hand,  were  affected  by  the  test  environment, 
thereby  reducing  delay  for  the  more  aggressive  environments. 

1.65  cm  (0.65  in.)  Thick  Plate 


To  further  examine  the  influences  of  environment  and  K level 
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on  fatigue  crack  growth  response  associated  with  a single  high- 
load  excursion  (overload) , experiments  were  carried  out  in  dehumid- 
ified argon,  air  and  3.5  pet  NaCl  solution  on  the  1.65  cm  thick 

plate.  Two  stress  intensity  levels  = 8.8  and  15.4  MN-m”^'^^ 

1/2 

or  8.0  and  14.0  ksi-in.  ' ),  with  an  overload  ratio  (K,  /K_  ) 

imax  2max 

of  2.0,  were  used.  Fatigue  crack  growth  was  carried  out  under 
sinusoidal  loading  at  a stress  ratio  of  0.05  and  a frequency  of  20 
Hz.  The  "constant  K"  procedure  was  utilized,  and  the  entire  experi- 
ment was  fully  automated.  Six  over-load  tests  were  run  on  each 
specimen.  Data  associated  with  each  over-load  test  are  designated 
by  a number  from  1 to  6;  the  numbers  correspond  to  increasing 

Experimental  data  for  tests  in  dehumidified  argon  are  shown 
in  Figure  11;  those  for  tests  in  air,  in  Figures  12  and  13;  and 
those  for  tests  in  3.5  pet  NaCl  solution  in  Figures  14  and  15. 
Because  of  more  precise  control  afforded  by  the  automated  testing 
procedure,  there  was  considerably  less  scatter  in  these  results  as 
compared  with  those  on  the  7.6  cm  (3.0  in.)  thick  plate  material 
(see  Figures  5 to  8 for  example) . The  results  show  that  the  over- 
load affected  zone  extended  out  a distance  of  about  2 r^  , from  the 
point  of  overload  application  and  was  essentially  independent  of 
the  test  environment.  The  data  also  indicate  that  the  delayed  re- 
tardation region  was  considerably  smaller  than  the  monotonic  plastic 
zone  size  order  of  the  cyclic  plastic  zone 

size  . Delay  in  the  aggressive  environments  (air  and  3.5 

pet  NaCl  solution)  was  shorter  than  in  dehvimidif ied  argon,  although 
there  was  no  discernible  difference  between  delay  (or  the  rate  of 
crack  growth)  in  air  and  in  3.5  pet  NaCl  solution  for  the  two 
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^2inax  Longer  delay  was  observed  at  the  lower 
^2max  both  environments  (Figures  13  and  15) . The  ob- 
served trend  on  delay  is  in  agreement  with  previous  results  on 
Ti-6A1-4V  alloy  [21-24]. 


DISCUSSION 


Data  Scatter 

Considerable  eunounts  of  scatter  were  observed  in  the  data 
for  specimens  taken  from  the  7.6  cm  (3.0  in.)  thick  plate,  par- 
ticularly for  those  obtained  under  "constant  X"  conditions  (see 
Figures  5 and  6,  for  example) . Several  sources  contributed  to 
this  scatter;  all  of  which  are  difficult  to  quantify  but  are  in- 
terrelated. First,  the  dc  electrical  potential  system  for  moni- 
toring crack  growth  readily  responds  to  external  disturbances, 
such  as  interruptions  in  load  and  changes  in  temperature,  and  has 
lower  noise  rejection  capabilities  than  the  ac  system.  Interrup- 
tions in  testing  to  apply  overloads  tended  to  perturb  the  system, 
which  required  some  time  (order  of  minutes)  to  re-establish  equi- 
librium. During  these  periods  of  perturbation,  the  potential  mea- 
surements can  be  in  error  by  several  tenths  of  microvolts  (corre- 
ponding  to  errors  in  crack  lengths  of  say  0.015  cm  or  0.006  in. 
for  0.3uV  disturbance  in  the  1.27  cm  thick  specimens,  and  lesser 
amounts  in  the  thinner  samples) . Errors  in  crack  growth  rates 
may  reach  the  order  of  ± 50  pet.  Second,  the  testing  procedure 
itself  can  be  sources  for  scatter.  The  manual  interruption  proce- 
dure used  in  conjunction  with  overload  application  and  the  resump- 
tion of  fatigue  required  "manual"  setting  and  adjustment  of  loads. 
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Each  load  change  required  time,  which  cannot  always  be  controlled. 
The  variations  in  these  "setup"  times  contributed  to  scatter  in 
two  ways;  (a)  It  altered  the  period  of  perturbation  for  the 
crack  monitoring  system,  and  thereby  contributed  to  scatter  as 
discussed  previously;  and  (b)  the  time  at  minimum  load  (and  possi- 
bly under  slight  compressive  load)  could  alter  the  subsequent 
crack  growth  response  [21,22].  Other  sources  of  variability  in- 
cluded accuracy  in  estimating  K during  testing  (specifically  for 
computing  and  at  the  overload  points)  , nonconstancy 

of  K,  the  influence  of  crack  length  on  the  rate  of  change  of  K 
with  crack  growth,  etc.  Overall,  uncertainty  in  Aa„  and  da/dN, 

Uij 

immediately  following  the  overload,  can  vary  by  the  cimounts  esti- 
mated previously. 

Because  of  the  use  of  automated  control  and  of  the  improved 
ability  for  rejecting  external  disturbances  with  the  ac  potential 
system,  scatter  in  data  from  the  1.65  cm  (0.65  in.)  thick  plate 
was  much  smaller.  Uncertainty  in  was  estimated  to  be  better 

than  0.005  cm  (0.002  in.)  for  the  1.65  cm  (0.65  in.)  thick  speci- 
mens. Uncertainty  in  crack  growth  rates  was  estimated  to  be  about 
± 20  pet.  Greater  reliance,  therefore,  will  be  placed  on  the  data 
from  the  1.65  cm  thick  plate  in  the  interpretation  of  experimental 
results . 

In  addition  to  these  sources  of  scatter,  the  proximity  of 
the  back  edge  of  the  specimen  at  long  crack  lengths  (greater  than, 
say,  4 cm)  may  have  also  affected  the  test  results.  This  possible 
effect  is  suggested  by  data  on  the  1.65  cm  thick  specimens  (see 
data  set  labeled  as  "6"  in  Figures  11  to  15) . Such  an  influence. 
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however,  was  not  apparent  in  the  data  on  specimens  machined 
from  the  7.6  cm  thick  plate  (see  Figures  5 and  6);  although  the 
effect  might  have  been  obscured  by  scatter  introduced  from  other 
sources.  Further  tests  would  be  needed  to  better  establish  the 
existence  and  significance  of  this  effect. 

Overload  Affected  Zone  and  Delayed  Retardation  Zone 

In  spite  of  the  scatter  in  some  of  the  data,  a reasonably 
consistent  trend  is  evident  in  all  of  the  data,  and  conforms  to 
that  illustrated  schematically  in  Figure  1.  The  high- load  excur- 
sion (overload)  affected  a zone  of  material  ahead  of  the  crack 
tip.  This  overload  affected  zone  is  delineated  by  the  point  of 
overload  application  to  the  point  where  steady-state  fatigue  crack 
growth  is  re-established  (that  is,  the  distance  from  Point  B to 
Point  E in  Figure  lb) . Within  the  overload  affected  zone,  fa- 
tigue crack  growth  response  included  (a)  an  abrupt  increase  in 
growth  rate  to  a value  greater  than  the  steady-state  rate  (or 
crack  acceleration) , (b)  a fairly  rapid  decrease  in  rate  to  a 

minimum  value  (delayed  retardation)  over  some  small  distance  ahead 
of  the  crack-tip  (delayed  retardation  zone) , and  (c)  subsequent 
gradual  increase  in  rate  back  to  the  steady-state  value. 

Examination  of  Figures  5,  7,  9a,  10a,  11a,  12  and  14  suggests 
that,  for  the  2219-T851  aluminum  alloy,  the  overload  affected  zone 
is  essentially  equal  to  the  monotonic  plastic  zone  size  associated 
with  the  overload  K,  that  is  see  Eqs.  (5)  and  (6)  [8,10,11]*. 

♦Because  of  the  contribution  of  surface  shear  lips  to  delay  [22, 

25],  the  lower  bound  for  Zql  be  larger  than  rjpi,  and  appears 
to  be  equal  to  about  2 rjpi  for  this  2219-T851  aluminum  alloy 
over  the  range  of  specimen  thickness  and  K^max  values  used  in 
this  investigation. 
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Depending  on  specimen  thickness  and  the  appropriate  zone 

size  may  range  between  that  for  "plane-strain"  and  that 

for  "plane-stress"  • The  zone  sizes  are  indicated  on  the 

various  figures.  The  overload  affected  zone  size  thus  obviously 
depends  on  and  on  specimen  thickness,  but  does  not  appear 

to  be  dependent  on  the  test  environment  (see  Figure  10  and  com- 
pare Figures  11a,  12  and  14) . Since  the  overload  affected  zone 
resulted  principally  from  mechanical  interactions,  the  apparent 
independence  of  zone  size  on  test  environment  is  reasonable. 

Significant  crack  acceleration  was  observed  in  some  cases 
and  not  in  others,  and  was  more  prevalent  for  tests  in  dehumid- 
ified argon  (an  inert  environment)  than  in  air  and  in  3.5  pet 
NaCl  solution  (aggressive  environments),  see  Figures  5 and  15. 

No  systematic  behavior  was  obvious,  and  the  crack  acceleration 
phenomenon  would  require  further  exeimination. 

Delayed  retardation  appeared  to  take  place  over  an  increment 
of  crack  growth  (Aa^^^)  that  ranged  from  the  cyclic  plastic  zone 
size,  Eqs.  (9)  and  (10) , to  the  monotonic  plastic  zone  size,  Eqs. 

(7)  and  (8),  associated  with  the  fatigue  loading  (that  is,  AK2 
or  . Depending  on  the  stress  intensity  level  and  specimen 

thickness,  again,  the  appropriate  zone  size  may  be  either  that  for 
"plane-strain"  or  that  for  "plane-stress".  In  light  of  the  smaller 
eimount  of  scatter  in  data  on  the  1.65  cm  thick  plate,  one  might 
suggest  that  the  delayed  retardation  zone  size  be  considered  to  be 
equal  to  the  cyclic  plastic  zone  size.  The  existence  of  delayed 
retardation  then  may  be  interpreted  in  terms  of  the  growth  of  the 
crack  through  a region  that  had  already  experienced  cyclic  fatigue 
damage. 


Effects  of  Specimen  Thickness,  Environment  and  K Level  on  Delay 

Having  established  that  the  overload  affected  zone  depended 
on  specimen  thickness  and  K level  (that  is,  on  crack-tip  con- 
straint) and  was  essentially  independent  of  test  environment,  the 
influences  of  these  variables  on  delay  (or  N*)  may  now  be  consid- 
ered. Delay  in  essence  represents  the  number  of  fatigue  load 
cycle  that  is  required  to  extend  the  crack  through  the  overload 
affected  zone.  As  such,  it  is  proportional  to  the  zone  size  and 
is  inversely  proportional  to  the  average  rate  of  crack  growth 
through  the  overload  affected  zone. 

For  specimens  of  different  thicknesses,  the  overload  affected 
zone  sizes  are  different,  with  larger  zone  sizes  associated  with 
the  thinner  specimens  (that  is,  those  tending  toward  "plane-stress") . 
For  the  range  of  K levels  used  in  investigation,  the  rate  of  fa- 
tigue crack  growth  in  a given  environment  was  sensibly  independ- 
ent of  specimen  thickness  [36,37].  Hence,  the  thinner  specimens 
would  tend  to  experience  greater  delay  (see  Figures  6 and  8) . 

This  observation  is  consistent  with  results  on  other  aluminxam 
alloys  [25,26]. 

A similar  argviment  can  be  applied  in  considering  the  effects 
of  environment.  Since  the  overload  affected  zone  size  does  not 
alter  with  environment,  delay  (N*)  would  simply  vary  inversely  with 
the  fatigue  crack  growth  rate.  Delay  in  the  more  aggressive  envi- 
ronment, therefore,  would  be  expected  to  be  shorter,  which  is  con- 
sistent with  observations  (see  Figures  10b,  11b,  13  and  15)*. 

*Note  that  in  Figures  13  and  15,  there  was  little  difference  in  the 
rate  of  fatigue  crack  growth  and  in  delay  between  tests  in  air  and 
in  3.5  pet  NaCl  at  a loading  frequency  of  20  Hz. 
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Interpretation  of  the  effects  of  K level  is  less  straight- 
forward because  increasing  K increases  both  the  overload  affected 

zone  size  and  the  rate  of  fatigue  crack  growth.  The  net  effect, 
however,  resulted  in  a decrease  in  delay  in  this  material  (see 
Figures  9b,  13  and  15) . This  observation  is  consistent  with  pre- 
vious results  on  a Ti-6A1-4V  alloy  [21-24],  but  differed  with  some 
results  on  a 7075-T6  aluminum  alloy  [25] . The  difference  presum- 
ably may  be  attributed  to  the  use  of  instead  of  N*  in  the  pre- 
vious experiments  [25]  (see  Figure  1) . Data  from  this  investiga- 
tion suggest  that  may  be  less  than  one-half  of  N*  when  there  is 
significant  crack  acceleration  following  an  overload. 


Modeling  Considerations 

A number  of  models  have  been  suggested  to  account  for  the 
effects  of  high-load  excursions  (overloads)  on  the  subsequent 
fatigue  crack  growth  response  [7-15].  A suitable  model,  however, 
should  include  a description  of  crack  acceleration  and  delayed 
retardation,  as  well  as  one  for  the  recovery  to  steady-state  crack 
growth  rate.  Since  the  modified  Willenborg  et  al.  model  [10-12] 
appeared  to  provide  useful  results,  a comparison  of  this  model 
with  the  present  results  seemed  appropriate.  Specifically,  this 
comparison  pertained  to  its  ability  to  predict  the  recovery  stage 
of  crack  growth  since  the  model  does  not  account  for  crack  accel- 
eration and  delayed  retardation.  Because  both  constant- amplitude 
fatigue  crack  growth  and  delay  data  are  availaUble  for  dehxjmidified 
argon,  a comparison  for  this  environment  was  made.  Specifically, 


data  obtained  on  a 1.65  cm  thick  specimen  at  K2max  * ^5.4 


MN-m 


-3/2 


(14  ksi-in 


1/2 


),  Figure  11,  were  used,  and  comparison  with  model 


predictions  is  shown  in  Figure  16. 
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According  to  the  modified  Willenborg  et  model  [10-12], 

the  stress  intensity  factors  for  fatigue  are  to  be  reduced  by  a 

residual  stress  intensity  factor  that  decayed  linearly  with 

crack  extension.  In  terms  of  notations  used  here,  K„  is  given 

£\ 

as  follows: 


M (1  - ) 

* T,W  _ th  2max 

S = ♦ ’'r  * 


Imax 


2max 


is  the  maximum  stress  intensity  factor  associated  with  fa- 
tigue crack  growth  threshold  at  R = 0;  Aa.  is  crack  growth  fol- 
lowing  the  overload;  and  is  the  overload  affected  zone* . S is 
defined  as  a shut-off  ratio,  and  corresponds  to  that  value  of  the 
ratio  of  K]iniax'^^2max  crack  arrest  is  expected  to  result. 

Taking  = 0.312  cm  (0.124  in.)  and  an  assumed  value  of 

K..  of  4.4  MN-m“3/2  (4.0  ksi-in.^/^)  [12,36-39],  values  of  K_  were 
■tn  K 

computed  for  two  shut-off  ratios  (that  is,  S = 2.3  and  2.8)**. 
These  values  of  K_  were  used  to  compute  effective  stress  intensity 

X\ 

factors  at  different  values  of  Aa^^^,  which  were  used  in  turn  to 
determine  the  corresponding  fatigue  crack  growth  rates  from  con- 
stant-amplitude data  [36-39].  The  model  predictions  are  shown  as 
dashed  lines  in  Figure  16. 

Inspection  of  Figure  16  indicates  that  better  agreement  was 
obtained  by  using  S = 2.8.  The  modified  Willenborg  et  a^.  model. 


* Zql  was  taken  to  be  (y/2it)  (K^jnax/'^ys)  » where  y is  an  empirical 
constant  [12]. 

**Shut  off  ratio  of  2.3  corresponded  to  the  value  used  by  Gallagher 
and  Stalnaker  [12],  and  S = 2.8  corresponded  to  data  on  a Ti- 
6A1-4V  alloy  [21,22] . 
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to  the  overload  plastic  zone  size  ('x-2 


however,  appeared  to  underestimate  the  rate  of  recovery  to  steady- 
state  growth.  This  observation  suggests  that  a faster  decay  in 

would  be  required.  Furthermore,  because  the  model  presumed  that 
the  effect  of  overload  was  exhausted  when  the  "leading  edge"  of  the 
fatigue-crack  plastic  zone  reached  the  boundary  of  the  overload 
plastic  zone,  an  overload-affected-zone  that  was  larger  than  the 
overload  plastic  zone  (in  this  case,  assumed 

to  model  the  recovery  response  (see  Figure  16).  The  results,  on 
the  other  hand,  suggested  that  the  overload  affected  zone  was  equal 

The  exhaustion  of  over- 
load effect  would  therefore  correspond  to  the  emergence  of  the  fatigue 
crack  from  the  overload  plastic  (or  overload  affected)  zone,  which 
provides  certain  physical  appeal. 

Based  on  these  observations  and  on  the  experimental  data,  it 
appeared  reasonable  to  suggest  further  semi-empirical  modifications 
to  the  model  of  Willenborg  et  a^.  [8] . In  these  further  modifi- 
cations, the  basic  concept  of  a residual  stress  intensity  factor 

produced  by  the  overload  is  preserved.  The  rate  of  decay  of 

2 

is  assumed  to  be  proportional  to  (1  - Aa^j^/Z^j^)  over  a range 
of  Aa^^  from  Z*^^  to  viz.: 


S ' ^OL  ^ ^ ^OL 


(12)' 


Z*j^  is  the  delayed-retardation-zone  size  and  is  assumed  to  be 
equal  to  the  appropriate  cyclic  plastic  zone  size  (or  plastic  zone 
size) . Z--  is  the  overload-affected-zone  size  and  is  taken  to  be 
equal  to  the  appropriate  overload  plastic  zone  size.  is  the 

residual  stress  intensity  factor  at  Aa„  * 0 (that  is,  immediately 
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following  the  overload)  and  is  given  by  [8,12]. 


- ‘'th/‘'2r,ax> 
(S  - 1) 


^^Imax  ^2max^ 


- ''th''''2max>T§-^  ''2max 


Q = K,  /K_  > 1 is  the  overload  ratio.  Crack  arrest  is  expect- 

Imax'^  2max 

ed  to  occur  when  Q i S [12] . 

According  to  this  modified  model,  a minimum  in  crack  growth 
rate  following  an  overload  would  be  expected  to  occur  at  = 

Z*  and  assume  a value  corresponding  to  the  effective  stress  in- 
tensity  factor  at  that  point,  as  determined  from  Eq.  (12)  . De- 
layed retardation  response  may  be,  as  a first  order  approximation, 
obtained  by  joining  a straight  line  between  the  point  of  overload 
application  and  the  minimum  growth  rate  point  (at  Aa^^^  = = 

r^p^) , and  could  be  expressed  in  exponential  form,  if  desired. 
Recovery  in  crack  growth  would  be  described  in  a manner  similar 
to  the  original  model  by  using  K_  given  by  Eq.  (12) . Estimated 
response  based  on  this  modified  model  (neglecting  crack  accelera- 
tion) is  indicated  by  the  solid  line  in  Figure  16  for  S = 2.8, 
and  is  seen  to  be  in  good  agreement  with  the  data.  Comparable 
agreement  was  obtained  with  data  from  the  7.6  cm  thick  plate  (that 


is,  data  shown  in  Figures  5 and  6) . Further  development  of  this 
empirical  model  and  experimental  verification  are  needed. 

SUMMARY 

To  expand  on  the  phenomenological  understanding  of  the 
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effects  of  load  interactions  on  fatigue  (principally  delay) , 
the  influences  of  plate  thickness,  K level  and  chemical  environ- 
ment on  fatigue  crack  growth  response  following  a single  high-load 
excursion  (overload)  were  examined  for  a 2219-T851  aluminum  alloy. 
An  overload  ratio  (Kj^inax'^^2max^  used.  The  principal 

results  of  this  investigation  are  as  follows: 

1.  The  high-load  excursion  (overload)  affected  a region 

of  material  ahead  of  the  crack  tip,  within  which  fatigue  crack 
growth  deviated  from  its  steady-state  behavior.  The  size  of  the 
overload  affected  zone  was  dependent  on  the  crack-tip  constraint, 
and  was  thereby  dependent  on  K level  and  specimen  thickness.  For 
a given  K level  and  specimen  thickness,  the  overload  affected  zone 
for  the  2219-T851  alloy  was  independent  of  chemical  environment  and 
was  found  to  be  about  equal  to  the  appropriate  (plane-strain  or 
plane-stress)  plastic  zone  size  (diameter)  for  the  high  load  (or 
^Imax^  (viz.,  2 for  "plane-strain"  and  r^^^  for  plane-stress"). 

2.  Delay  (N*)  was  found  to  increase  with  decreasing  specimen 
thickness.  This  increase  in  delay  may  be  attributed  primarily  to 
an  enlargement  in  the  overload-affected-zone  size  for  the  thinner 
specimens . 

3.  Delay  (N*)  was  found  to  decrease  with  increasing  aggres- 
siveness of  the  environment  (that  is,  dehumidified  argon  air  -*• 

3.5  pet  NaCl  solution).  This  decrease  in  delay  may  be  attributed 
to  the  higher  rates  of  fatigue  crack  growth  in  the  more  aggressive 
environments . 

4.  Decreasing  K level  produced  an  increase  in  delay  in  this 
alloy.  This  dependence  on  K was  the  net  result  of  the  separate 
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influences  of  K on  the  overload-affected-zone  size  and  on  the 
rate  of  fatigue  crack  growth. 

5.  The  results  provided  further  confirmation  of  delayed 
retardation,  and  suggested  that  this  phenomenon  is  associated  with 
crack  growth  through  a region  of  material  within  the  overload  af- 
fected zone  that  had  sustained  prior  fatigue  dcimage.  Based  on  this 
interpretation,  it  is  suggested  that  a delayed  retardation  zone  be 
defined,  and  that  the  size  of  this  zone  be  identified  with  the 
cyclic  plastic  zone  size  associated  with  steady-state  fatigue  im- 
mediately preceding  the  overload.  Minimum  in  the  rate  of  fatigue 
crack  growth  within  the  overload  affected  zone  may  be  expected  to 
occur  at  the  delayed  retardation  zone  boundary.  These  suggestions 
tended  to  be  supported  by  the  data,  although  some  of  the  data  sug- 
gested that  this  zone  size  may  be  as  large  as  the  monotonic  plas- 
tic zone  diameter  associated  with  prior  fatigue. 

6.  A residual  stress  intensity  factor  (Kj^)  concept  appeared 
to  be  useful  in  describing  fatigue  crack  growth  response  following 
an  overload  application.  A further  modification  to  the  model  of 
Willenborg  et  al . , to  increase  the  rate  of  decay  in  was  sug- 
gested. This  modification  appeared  to  provide  reasonable  estimates 
of  fatigue  crack  growth  response. 
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Table  1 


Chemical  Composition  and  Tensile  Properties  of  7.6  cm  (3-in.) 
thick  2219-T851  Aluminum  Alloy  Plate* 

a.  Chemical  Composition  (Weight  Percent) 

Cu  Zn  Mg  Si*  Mn  Ti  Cr  Fe  A1 
6.28  0.025  0.003  0.088  0.25  0.051  <0.0001  0.25  Balance 


Orientation 


b.  Tensile  Properties** 

Yield  Tensile  Elongation 

Strength  Strength  in  5.08  cm 

MN/m^ (ksi)  MN/m^ (ksi)  pet 


Reduction 
of  Area 


Longitudinal  358  (52) 


455  (66) 


Transverse 


351  (51) 


455  (66) 


* Data  from  Westinghouse  Electric  Corporation. 

**  Averages  of  6 tests  from  2 locations. 


Notes;  All  dimensions  in  cm  ( 2.54cm  = I ia) 

C 3 current  leads',  potential  leads 

Figure  2:  Wedge-opening- load  (WOL)  specimens  and  place- 
ment of  current  and  potential  leads  for  the 
crack  monitoring  system. 


(corrected)  (in.) 


NORMALIZED  POTENTIAL.  V* • (V(OI-Vr)/Vr 


(a) 


normalized  potential,  V**(V(0l-Vr)/V, 


(b) 

Figure  3;  Calibration  curve  for  crack  length  versus 
normalized  dc  potential  for  specimens  from 
(a)  the  7.6  cm  thick  plate  (see  Figure  2a), 
and  (b)  the  1.65  cm  thick  plate  (see  Figure 
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Aool*  0 - Oql  (in  ) 

002  0 04  006  008  0 10  012 


-012  -008  004  0 004  008  0.12  016  0 20  024  0 28 


Aa^t*  0 - OoL  (cm) 

(a) 


Aag^>  0 - Og^  (in  ) 


AOg^"  0 - Og^  (cm) 


(b) 


Figure  5:  Fatigue  crack  growth  response  (da/dN  versus 
AaoL  following  a single  overload  at  Kimax/ 
K2max  = 2.0  using  the  "constant  K"  procedure; 
(a)  0.25  cm  thick  specimens,  and  (b)  1.27  cm 
thick  specimens.  (Test  environment:  dehumid- 
ified argon  at  room  temperature.  Fatigue 
loading  at  R » 0.05). 


•>/«N  |un/cyclt)  *’/dN  l“'/e»ct«) 


-2  0 2 4 6 8 10 

ANot«N-  Not(IO«ercl«) 


(a) 


(b) 

Figure  6:  Fatigue  crack  growth  response  (da/dN  versus 
ANol)  following  a single  overload  at  Klmax/ 
K2max  = 2.0  using  the  "constant  K"  procedure 
(a)  0.25  cm  thick  specimens,  and  (b)  1.27 
cm  thick  specimens.  (Test  environment: 
dehumidified  argon  at  room  temperature. 
Fatigue  loading  at  R = 0.05). 
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Figure  7;  Fatigue  crack  growth  response  (da/dN  versus 
AaoL)  following  a single  overload  at  Kimax/ 
K2max  = 2.0  using  the  "constant  load"  proce- 
dure; (a)  0.51  cm  thick  specimen,  and  (b) 
1.27  cm  thick  specimen.  (Test  environment; 
dehumidified  argon  at  room  temperature. 
Fatigue  loading  at  R **  0.05). 
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Figure  8:  Fatigue  crack  growth  response  (da/dN  versus 
ANql)  following  a single  overload  at  Kimax/ 
K2max  “ 2.0  using  the  "constant  load"  proce 
dure:  (a)  0.51  cm  thick  specimen,  and  (b) 
1.27  cm  thick  specimen.  (Test  environment: 
dehumidified  argon  at  room  temperature. 
Fatigue  loading  at  R * 0.05). 


Figure  9;  Fatigue  crack  growth  response  for  a 0.51  cm 
thick  specimen  following  a single  overload 
^lmax/^2max  = 2.0  using  the  "constant 
load*  procedure:  (a)  da/dN  versus  AaoL»  and 
(b)  da/dN  versus  ANql*  (Test  environment: 
air  at  room  temperature.  Fatigue  loading  at 
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Figure  10:  The  influence  of  chemical  environment  on 
the  room  temperature  fatigue  crack  growth 
response  following  a single  overload  at 
Klmax/K2max  * 2.0,  with  K2inax  • 10.8 
MN-m~3/2  and  R = 0.05:  (a)  da/dN  versus 
AaoL»  and  (b)  da/dN  versus  ANql* 


dO/UN  ('"’/cycle)  *’/((N  ('"'/cycle) 
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(b) 

Figure  11;  Fatigue  crack  growth  response  for  a 1.65  cm 

thick  specimen  following  a single  overloaci  at 
^lmax/^2max  “ 2.0  and  K2max  = 15.4  MN-m~3/2 
using  the  "constant  K"  procedure:  (a)  da/dN 
versus  AaoL»  da/dN  versus  ANql. 

(Test  environment;  dehumidified  argon  at 
room  temperature.  Fatigue  loading  at  R = 
0.05)  . 
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Figure  12;  Fatigue  crack  growth  response  (da/dN  versus 
Aa©^)  for  a 1.65  cm  thick  specimen  following 
a single  overload  at  Kiniax/K2max  = 2.0  using 
the  "constant  K"  procedure:  (a)  K2max  = 8-8 
MN-m-3/2,  and  (b)  K2max  = 15.4  MN-m-3/2. 
(Test  environment;  air  at  room  temperature. 
Fatigue  loading  at  R * 0.05). 
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Figure  13:  Fatigue  crack  growth  response  (da/dN  versus 
ANqi,)  for  a 1.65  cm  thick  specimen  following 
a single  overload  at  Kimax/K2max  = 2.0  using 
the  "constant  K"  procedure':  (a)  K2max  = 5.8 
MN-m"3/2,  and  (b)  K2max  = 15.4  MN-m-3/2. 
(Test  environment:  air  at  room  temperature. 
Fatigue  loading  at  R = 0.05). 
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Figure  14:  Fatigue  crack  growth  response  (da/dN  versus 

for  a 1.65  cm  thick  specimen  following 
a single  overload  at  Kin,ax/K2max  = 2.0  using 
the  "constant  K"  procedure:  (a)  K2inax  = 8-8 
MN-m-3/2,  and  (b)  K2inax  = 15.4  MN-m-3/2. 
(Test  environment;  3.5  pet  NaCl  solution  at 
room  temperature.  Fatigue  loading  at  R = 
0.05). 
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Figure  15:  Fatigue  crack  growth  response  (da/dN  versus 
ANql)  for  a 1.65  cm  thick  specimen  following 
a single  overload  at  Kiinax/^2max  ” using 
the  "constant  K"  procedure:  (a)  K2inax  ® 8.8 
MN-m“3/2,  and  (b)  K2max  = 15.4  MN-m-3/2. 
(Test  environment;  3.5  pet  NaCl  solution  at 
room  temperature.  Fatigue  loading  at  R * 
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Figure  16;  Comparison  of  fatigue  crack  growth  response  with  "predictions"  based 
on  residual  stress  intensity  factor  concept. 
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CONSIDEi^^TIONS  OF  MEASUREMENT  PRECISION  AND  DATA  PROCESSING 
PROCEDv  ;S  FOR  FATIGUE  CRACK  GROWTH  TESTING  AND  ANALYSIS 

by 

R.  P.  Wei  and  W.  Wei* 

LEHIGH  UNIVERSITY 
Bethlehem,  Pennsylvania  18015 

ABSTRACT 

A computer  simulation  of  the  primary  fatigue  crack  growth 
data  (crack  length  versus  elapsed  cycles  or  a versus  N)  was  made 
to  examine  the  influences  of  crack  length  measurement  interval 
and  precision,  and  of  the  associated  data  processing  procedure 
on  the  resulting  data  on  fatigue  crack  growth  kinetics  (i.e., 
da/dN  versus  AK) . Variability  in  the  derived  growth  rate  (da/dN) 
data  depended  strongly  on  the  magnitude  of  the  measurement  inter- 
val relative  to  the  measurement  precision.  It  was  reduced  by 
those  data  processing  procedures  (such  as  the  incremental  poly- 
nomial methods)  that  fitted  a smooth  curve  through  portions  of  the 
primary  data.  Such  procedures,  however,  introduced  significant 
bias  into  the  derived  data  on  fatigue  crack  growth  kinetics.  The 
results  further  suggested  that  much  of  the  variability  in  the  pub- 
lished data  on  fatigue  crack  growth  kinetics  might  be  attributed 
to  sources  described  herein.  Caution  is,  therefore,  recommended 
in  attempting  to  draw  statistical  inferences  regarding  material 
variability  from  these  data. 


*Now  Research  Assistant  in  the  Department  of  Metallurgy  and  Mining 
Engineering,  University  of  Illinois  at  Urbana-Champaign,  Urbana, 
Illinois  61801. 


INTRODUCTION 


Many  engineering  structures  are  subjected  to  cyclically 
varying  (fatigue)  loads  during  service.  Fatigue  analysis,  or 
the  estimation  of  fatigue  lives  of  such  structures,  is  there- 
fore an  essential  part  of  engineering  design.  Fracture  mech- 
anics technology  [1-5],  developed  over  the  past  twenty  years,  is 
well  established  as  an  important  tool  for  fatigue  crack  growth 
analysis  and  for  the  analysis  of  related  structural  reliability 
problems.  Standardized  test  methods  have  been,  or  are  being, 
developed  to  obtain  supporting  experimental  data  for  these  anal- 
yses (for  example,  [6,7]).  For  fatigue  crack  growth,  a typical 
design  sequence  would  involve  the  following  steps  [5]  which  are 
illustrated  schematically  in  Figure  1: 

1.  The  acquisition  of  fatigue  crack  growth  data  (crack 
length  versus  number  of  elapsed  cycles,  or  a vs.  N) 
from  laboratory  specimens  or  from  model  components. 

2.  Processing  these  data  to  convert  them  into  fatigue 
crack  growth  rate  (da/dN  or  Aa/AN)  as  a function  of 
the  mechanical  driving  force.  The  mechanical  driving 
force  is  characterized  by  the  fracture  mechanics  para- 
meters AK  or  K (where  AK  and  K are  the  stress 

indx  ludx 

intensity  factor  range  and  maximum  stress  intensity 
factor,  respectively,  for  a given  load  cycle  [1,2], 
and  are  functions  of  the  applied  load  and  current 
crack  length  [5,7]. 

3.  Modeling  or  expressing  the  fatigue  crack  growth  ki- 
netics (i.e.,  the  da/dN  vs.  AK  or  K data)  by  an 

luciX 


r 
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appropriate  analytical  or  empirical  relationship. 

4.  Utilizing  the  analytical  or  empirical  relationship 
thus  derived  to  predict  or  estimate  the  fatigue 
crack  growth  response  for  the  structure  and  the  ex- 
pected service  conditions  of  interest. 

As  a part  of  the  prediction  procedure,  an  estimate  of  struc- 
tural reliability  would  be  required.  This  estimate  is  provided 
through  an  estimate  of  the  probable  variations  in  fatigue  lives 
or  the  probability  for  survival  following  a given  period  of  ser- 
vice. Thus,  it  is  important  to  incorporate  into  the  analyses  the 
variability  in  fatigue  crack  growth  rate  data  and  estimates  of 
the  accuracy  of  various  analysis  steps.  Variability  in  fatigue 
data  arises  from  several  sources:  (a)  inherent  variations  in 
material  properties  (particularly  its  fatigue  crack  growth  re- 
sistance) , (b)  variations  in  loading  and  environmental  conditions 
during  testing,  and  (c)  variations  introduced  by  crack  length 
measurement  precision  and  the  associated  data  processing  proce- 
dures. Unfortunately,  however,  variability  introduced  by  these 
several  sources  are  not  clearly  recognized  and  reported  separately 
along  with  the  kinetic  data  on  fatigue  crack  growth.  As  a result, 
it  is  commonly  (though  incorrectly)  assumed  that  the  observed 
scatter  in  the  kinetic  data  resulted  principally  from  material 
property  variations,  and  an  upper-bound  curve  to  these  data  is  of- 
ten used  for  design  analyses  [4].  This  assumption,  however,  can- 
not be  fully  justified,  because  the  fatigue  data  are  often  ob- 
tained from  a small  number  of  specimens;  often  taken  from  a single 
plate  of  material.  Hence,  the  data  can  reflect,  at  best,  local 
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variations  in  material  properties,  and  do  not  contain  the  type 
of  statistical  information  (in  terms  of  heat-to-heat  or  plate-to- 
plate  variations)  that  would  be  required  for  design. 

In  this  investigation,  variability  in  fatigue  crack  growth 
rate  data,  introduced  by  crack  length  measurement  interval  and 
precision  and  by  the  associated  data  processing  procedures,  was 
examined  systematically  through  computer  simulation.  This  inves- 
tigation expanded  upon  preliminary  studies  by  Clark  and  Hudak  [4] 
and  by  Graham  [8]  on  the  effects  of  data  processing  procedure  on 
variability.  The  possible  influence  of  data  processing  procedures 
on  the  accurate  representation  of  fatigue  crack  growth  kinetics 
(i.e.,  in  introducing  bias)  was  also  exeunined.  Information  from 
this  investigation,  in  conjunction  with  assessments  of  the  influ- 
ences of  loading  and  environmental  variables  through  the  usual 
error  analysis  procedures  [9] , would  provide  a basis  for  deducing 
variability  introduced  by  material  property  variations  (the  re- 
quisite quantity)  from  the  experimental  fatigue  data. 

prcx:edure 

Because  the  primary  experimental  fatigue  crack  growth  (a  vs. 
N)  data  would  contain  contributions  from  measurement  precision,  as 
well  as  the  effects  of  possible  variations  in  material  and  of 
loading  and  environmental  variables,  it  was  not  practicable  to 
utilize  such  data  to  assess  the  influences  of  measurement  preci- 
sion and  of  the  associated  data  processing  procedures.  If  it  were 
possible  to  exclude  all  other  effects,  the  assessment  procedure 
would  have  involved  reconstruction  of  the  fatigue  crack  growth 
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(a  vs.  N)  curve  {by  following  steps  1 through  4 described  pre- 
viously and  illustrated  in  Figure  1)  and  comparison  of  this 
curve  with  the  original  crack  growth  data.  Since  this  was  not 
possible,  an  alternative  procedure,  first  suggested  by  Clark 
and  Hudak  [4] , was  used.  In  this  procedure,  an  analytical  re- 
lationship was  assumed  for  the  fatigue  crack  growth  kinetics. 

For  simplicity,  a simple  exponential  relationship  was  assumed  to 
apply  over  the  entire  range  of  crack  lengths  of  interest  [4,10]. 

da/dN  = C(AK)"  (1) 

C and  n are  constants  (normally  to  be  determined  experimenta' 
ly) An  "error  free"  fatigue  crack  growth  (a  vs.  N)  curve  was 
then  created  from  this  kinetic  relationship  by  integrating  Equa- 
tion (1)  between  selected  initial  and  final  crack  lengths.  From 
this  curve,  an  "error  free"  a vs.  N data  set  is  produced  based  on 
the  chosen  crack  length  measurement  interval  (Aa) . Appropriate 
random  measurement  errors  were  then  generated  by  computer  simula- 
tion (for  a prescribed  measurement  precision)  and  were  applied  to 
the  "error  free"  a vs.  N data  set  to  produce  a set  of  simulated 
fatigue  crack  growth  (a  vs.  N)  data.  The  simulated  primary  (a  vs. 
N)  data  were  then  processed  by  the  various  data  processing  pro- 
cedures to  determine  the  fatigue  crack  growth  kinetics.  The  re- 
sulting kinetic  (derived)  data  were  finally  used  to  assess  vari- 
ability introduced  by  the  choice  of  measurement  intervals  and 

i/Note  that  C is  a dimensioned  constant,  whose  value  and  di- 
mension change  with  n to  provide  growth  rate  in  terms  of 
length  per  cycle. 
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measurement  precision.—'  Least  squares  error  procedure  was 
also  used  to  fit  the  derived  kinetic  data  to  an  equation  of  the 
form  given  by  Equation  (1) . The  constants  C and  n determined 
from  this  procedure  were  then  compared  with  the  original  C and  n 
values  to  determine  how  well  these  values  were  reconstructed,  and 
also  compared  amongst  the  different  methods  to  test  for  possible 
bias.  In  essence,  this  alternate  procedure  used  step  3 (illus- 
trated in  Figure  1)  as  a starting  point,  and  then  proceeded  through 
steps  4,  1 and  2 to  return  to  step  3.  By  this  simulation  proce- 
dure, the  effects  of  measurement  precision  were  isolated  for  study. 

A compact  specimen  [7],  having  a width  of  2.0  in.,  was 

used  as  a model  for  this  investigation.  The  constants  C and  n 

were  assumed  to  be  C = 4 x 10  ^ (in. /cycle)  (ksi-in}'*^^)  ” and  n = 

2.25  (see  Equation  (1))  [4].  Simulated  crack  growth  ranged  from 

0.6  to  1.6  in.  Stress  intensity  factor  solution  and  other  de- 
l/2 

tails  for  this  specimen  are  given  in  [7] . (Note  that  1 ksi-in. 

= 1.1  MN-m  and  1 in.  = 2.54  cm.) 

Unlike  Clark  and  Hudak  [4] , who  made  a simplifying  assump- 
tion that  the  measurement  error  was  constant  in  magnitude  and 
varied  alternately  in  sign  between  successive  measurement  points, 
measurement  error  was  assumed  to  vary  randomly  both  in  magnitude 
and  in  sign  in  this  investigation.  The  measurement  error  was 
assumed  to  be  described  by  a normal  distribution  function  with  a 


^Note  that  measurement  precision  would  affect  both  da/dN  and  AK. 
The  values  of  AK  would  be  altered  by  less  than  1 pet,  however, 
in  a typical  experiment,  and  need  not  be  considered  as  statisti- 
cally varying  quantities  in  the  analyses. 
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mean  of  zero  and  a known  (prescribed)  standard  deviation,  a. 
Measurement  error  for  each  successive  measurement  point  was  se- 
lected from  this  distribution  with  the  aid  of  a random  number 
generator.  Its  magnitude  was  rounded  to  a value  nearest  to  the 
assumed  resolution  of  the  measurement  instrument;  viz.,  to  0.40.^^ 
The  error  assignment  procedure  is  illustrated  in  Figure  2.  By 
this  procedure,  a more  realistic  simulation  of  the  random  experi- 
mental errors  was  achieved. 

The  following  measurement  intervals  and  measurement  preci- 
sion (given  in  terms  of  a)  were  selected  for  study.  The  measure- 
ment intervals  are  typical  of  those  in  common  use  and  are  given 
in  English  units  here.  (Note  that  1 in.  = 2.54  cm). 


Interval 

&a 

0.050  in. 
0.025  in. 
0.010  in. 


Precision 

a 

0.010,  0.005  and  0.0025  in. 
0.005,  0.0025  and  0.00125  in. 
0.002  and  0.001  in. 


For  each  combination  of  crack  length  measurement  interval  and 
precision,  ten  different  simulated  fatigue  crack  growth  (a,N) 
data  sets  were  generated  to  provide  sufficient  information  for 
statistical  analysis.  The  different  data  sets  were  produced  by 
using  randomly  selected  numbers  as  "seeds"  for  the  random  number 
generator,  and  simulated,  for  example,  data  that  might  be  obtain- 
ed from  testing  10  ideally  identical  specimens  by  a single 


When  an  infrequent  error  of  greater  than  2a  was  indicated,  a 
value  equal  to  2a  was  assigned. 
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experimenter  or  observations  of  10  different  experimenters  for  a 
single  specimen. 

For  this  investigation,  the  variability  factors  VF  and  VF* , 
defined  by  Clark  and  Hadak  [4,6],  were  adopted  for  measuring  vari- 
ability in  the  data  on  fatigue  crack  growth  kinetics.  VF  provided 
a measure  of  data  scatter  about  a least-squares-error  curve  through 
the  derived  data  (of  the  form  given  by  Equation  (1);  that  is, 
using  the  derived  values  for  C and  n) , whereas  VF*  provided  a • 
measure  of  scatter  against  the  original  curve.  VF  is  defined  by 
the  following  relationship  [4,6]. 

VF  = e^^  (2) 

R is  the  residual  standard  error  defined  as  follows: 


R = 


I [iln  (da/dN)  ^ - )ln(da/dN)^] 


Z - 2 


1/2 


(3) 


where 

(da/dN)j^  = "measured"  crack  growth  rate  at  AK^; 
(dS/dN)^  = "predicted"  crack  growth  rate  at  AK^, 
based  on  the  least-squares-error  fit; 

Z = number  of  data  pairs;  that  is  (da/dN)^ 
and  AK^  pairs. 

VF*  is  similarly  defined,  except  (dS/dN)^  would  be  given  by  the 
original  constants  C and  n. 

To  excunine  the  influences  of  crack  length  measurement  inter- 
val and  precision,  and  of  the  associated  data  processing  proce- 
dures on  variability,  two  data  processing  procedures  that  are  in 
common  usage  were  used  [6,7].  The  first  procedure  is  the  finite 
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difference  (or  secant)  method.  The  crack  growth  rate  was  ob- 
tained by  simply  taking  the  growth  increment  between  two  adjacent 
points  and  dividing  it  by  the  number  of  elapsed  cycles. 


O’*/®’!  = n' 

1+1  1 


(4) 


The  average  crack  length,  + a^)/2,  was  used  to  compute  the 

corresponding  value  for  AK^.  The  second  method,  the  seven-point 
incremental  polynomial  method,  involved  least-squares-error  fit- 
ting of  a second-order  polynomial  in  N to  seven  successive  (a,N) 
data  pairs  (i.e.,  over  the  range  < a < a^^^)  • Crack  growth 

rate  was  then  computed  as  a derivative  to  this  polynomial,  along 
with  AK^,  at  the  central  (or  ith)  point  [6,7].  The  following 
relations  were  used  in  calculating  da/dN  [6,7]; 

2 


a = b^  + b. 
o 1 


fN-cJ 

^2 

^ ^2 

""2 

(5) 


b (N  - C ) 

(da/dN)  . = + 2b,  ^ 

^ ^ ci 


(6) 


where 


C,  = 5 (N..3  . 


=2  - I '"ls  - "i-3' 


The  numbers  of  elapsed  cycles  were  normalized  in  terms  of  and 
C2  to  circumvent  the  need  of  dealing  with  large  numbers  in  the 
least-squares-error  fitting  routine  [4,6]. 

For  exaunining  the  influences  of  data  processing  procedures, 
the  incremental  polynomial  procedure  was  expanded  to  include  5 and 
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9 points.  In  addition,  a procedure  that  involved  fitting  para- 
bolas through  3 adjacent  data  points  (Simpson's  rule)  was  also 
included;  with  the  computation  of  growth  rate  made  again  at  the 
central  point.  A measurement  interval  of  0.025  in.,  with  a * 

0.005  and  0.00125  in.,  was  used  in  this  part  of  the  study.  Each 
data  processing  procedure  was  used  to  analyze  the  same  10  sets  of 
simulated  (a,N)  data  corresponding  to  the  chosen  measurement  in- 

1 

terval  and  precision.  Adjustments  were  made  to  ensure  that  the  i 

t 

same  number  of  (da/dN,AK)  data  pairs  (covering  the  same  range  of  i 

crack  lengths  or  AK)  was  included  in  each  procedure,  and  were 
achieved  by  excluding  appropriate  numbers  of  (da/dN,AK)  pairs  at 
the  beginning  and  at  the  end  to  yield  the  Scune  number  of  data 
pairs  as  that  of  the  9-point  incremental  polynomial  method. 

Analyses  of  the  results  were  made  by  comparing  the  mean 
values  of  C and  n from  each  data  processing  procedure (with  10 
sets  for  each  condition)  at  the  95  percent  confidence  level  [11, 

12] . The  statistical  tests  used  required  that  each  group  of  ten 
C and  n values  from  each  procedure  and  condition  be  normally  dis- 
tributed. All  of  the  C and  n values  were  assessed  [11,12]  and 
were  found  to  satisfactorily  fulfill  the  normality  requirement  in 
each  case. 

RESULTS  AND  DISCUSSIONS 

The  influences  of  measurement  interval  and  precision  and  of 
the  associate  data  processing  procedures  on  variability,  and  the 
possible  influences  of  the  various  data  processing  procedures  in 
introducing  bias,  are  reported  and  discussed  separately. 
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Measurement  Interval  and  Precision 

The  results  obtained  from  an  investigation  of  the  influences 
of  crack  length  measurement  interval  and  precision  and  of  the 
associated  data  processing  procedure  on  variability  in  fatigue 
crack  growth  rates  are  summarized  in  Table  1.  It  can  be  seen 
that  the  variability  (as  measured  by  VF  and  VF*)  decreased  as  the 
measurement  precision  was  increased  (i.e.,  for  decreasing  values 
of  a) . For  a given  measurement  precision  (or  a) , the  results  on 
variability  appear  to  be  somewhat  ambiguous.  For  the  finite  dif- 
ference (secant)  method,  the  variability  increased  with  decreasing 
measurement  interval.  For  the  incremental  polynomial  method, 
variability  (as  indicated  by  VF)  increased  with  decreasing  meas- 
urement interval,  whereas  VF*  appeared  to  indicate  improvement 
(from  Aa  = 0.050  in.  to  Aa  = 0.025  in.)  before  it  worsened.  For 
a given  ratio  of  a to  Aa  (for  excimple,  at  a/Aa  = 0.1)  there  ap- 
peared to  be  an  initial  slight  improvement  followed  by  little 
change  in  variability  with  decrease  in  measurement  interval.  One 
must  recognize  that  the  total  number  of  (a,N) , and  the  corre- 
sponding (da/dN,AK)  data  pairs  increased  with  decreasing  measure- 
ment interval;  the  number  of  (a,N)  data  points  being  20,  40  and 
100  respectively  for  Aa  of  0.050,  0.025  and  0.010  in.  in  the  pres- 
ent case.  Since  the  variability  is  associated  with  the  same 
population,  the  observed  trend  is  to  be  expected. 

Comparison  between  VF  and  VF*  for  the  two  data  processing 
procedures  indicated  that  the  secant  method  produced  greater  var- 
iability. This  result  was  expected,  because  the  incremental  poly- 
nomial method  tended  to  "smooth  out"  the  data.  Comparison 


between  VF  and  VF*,  for  a given  data  processing  procedure,  indi- 
cated that  VF  and  VF*  for  a given  La.,  a combination  were  essen- 
tially equal  for  the  secant  method,  whereas  they  were  significant- 
ly different  for  the  incremental  polynomial  method.  The  larger 
VF*  values  suggested  that  data  did  not  agree  with  the  originally 
assumed  curve,  and  that  the  incremental  polynomial  method  tended 
to  introduce  a systematic  bias  into  the  processed  data.  This  ob- 
servation is  supported  by  comparing  the  tabulated  values  of  C and 
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n against  the  originally  assumed  values  of  C *=  4 x 10  and  n * 
2.25.  A more  detailed  evaluation  of  the  data  processing  procedure 
is  given  in  the  following  section. 

Examination  of  VF  and  VF*  in  Table  1 (particularly  for  the 
secant  method)  shows  that  significant  variability  can  be  intro- 
duced by  the  injudicious  choice  of  measurement  interval  in  rela- 
tion to  the  measurement  precision.  It  appears  that  a measurement 
interval  equal  to  10  times  the  measurement  precision  (i.e. , La  = 

10  a)  would  produce  a "tolerable"  scatter  in  the  fatigue  crack 
growth  rate  data.  There  is,  however,  an  upper  limit  to  the  choice 
of  measurement  interval.  This  choice  is  dictated  by  the  desired 
precision  in  identifying  crack  length  and  the  corresponding  AK 
with  the  measured  crack  growth  rate.  The  choice  of  measurement 
interval  must  be  based  also  on  a realistic  assessment  of  measure- 
ment precision  (which  must  be  established  experimentally  and  is 
not  the  same  as  the  smallest  division  on  a measuring  instrument  — 
viz.  instrumental  resolution  — as  is  often  assumed) . Assessment 
of  measurement  precision  may  be  made,  for  exeunple,  from  independ- 
ent measurements  of  several  crack  specimens  by  different  investi- 


12 


gators,  or  from  repeated  measurements  of  these  same  specimens 
by  the  same  investigator,  taken  at  random  and  at  different  times 
to  minimize  bias. 

The  magnitudes  of  the  variability  factors  (VF  and  VF*)  given 
in  Tables  1 and  2 are  comparable  to  those  deduced  from  typical 
experimental  data  (see  [6],  for  example).  This  comparison  sug- 
gests that  much  of  the  scatter  in  the  published  data  on  fatigue 
crack  growth  kinetics  may  be  attributed  to  the  measurement  and 
data  processing  procedures.  Until  variability  from  these  sources, 
along  with  that  introduced  by  load  and  environment  variation  dur- 
ing each  test  and  from  test  to  test,  can  be  separated  from  the 
experimental  data  (obtained  from  sufficient  numbers  of  specimens) , 
no  meaningful  information  with  respect  to  material  property  var- 
iations can  be  deduced.  One  should,  therefore,  refrain  from 
drawing  statistical  inferences  with  respect  to  material  property 
variations  directly  from  the  experimental  data  on  fatigue  crack 
growth  kinetics  at  this  time. 

Data  Processing  Procedure 

C,  n,  VF  and  VF*  values,  along  with  the  estimates  of  stand- 
ard errors,  are  summarized  in  Table  2 and  are  given  in  two  parts. 
In  the  first  part,  the  values  were  deduced  by  using  the  seune 
number  of  points  (32),  over  the  scime  range  of  crack  lengths  or  AK, 
for  each  procedure.  In  the  second  part,  the  values  were  deduced 
from  all  of  the  data  points  produced  by  each  procedure.  Informa- 
tion from  the  first  part  was  used  to  compare  the  data  processing 
procedures  against  one  another,  and  the  results  are  siimmarized  in 
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Table  3.  The  second  part  was  used  to  compare  the  individual 
procedures  against  the  original  curve.  This  comparison  is  given 
in  Table  4. 

Statistical  tests  [11,12],  at  the  95  percent  confidence 
level,  showed  that  there  was  a systematic  bias  introduced  by  the 
incremental  polynomial  methods,  as  indicated  by  reduced  values  of 
C and  increased  n (see  Tables  2 to  4) . The  bias  increased  with 
increasing  number  of  points  used  in  the  incremental  polynomial 
method.  The  finite  difference  (secant)  method  and  the  parabolic 
curve  fitting  procedure  (Simpson's  rule)  did  not  show  significant 
bias,  although  these  methods  produced  greater  scatter  (see  Tables 
2 to  4) . It  appears  that  one  of  these  two  methods  should  be  used 
for  data  processing,  at  least  in  a fundamental  sense  (for  example, 
to  try  to  deduce  a fundamental  fatigue  crack  growth  relationship 
from  experimental  data).  In  a practical  sense,  the  bias  may  not 
be  significant  in  terms  of  probable  errors  in  predicted  fatigue 
lives.  The  question  of  practical  significance  needs  to  be  estab- 
lished by  more  detailed  studies  covering  a broad  range  of  da/dN 
versus  AK  relationships  and  a wider  range  of  crack  lengths. 

CONCLUSION 

The  results  of  this  investigation  indicate  that  considerable 
variability  in  fatigue  crack  growth  rate  data  can  result  from  the 
choice  of  crack  length  measurement  interval  in  relation  to  the 
measurement  precision,  and  the  associated  data  processing  proce- 
dures. Data  processing  procedures  that  tend  to  reduce  scatter, 
through  the  use  of  "smoothing"  techniques,  could  introduce  bias 
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into  the  results.  With  suitable  choice  of  measurement  interval 
in  relation  to  precision  (viz.,  Aa/o  = 10  to  20),  data  obtained 
by  the  secant  method  appears  to  best  represent  actual  response, 
with  an  acceptable  degree  of  variability  (i.e.,  VF  < 2) . One 
must  be  realistic,  however,  in  assessing  the  measurement  preci- 
sion. 

These  results  further  suggested  that  much  of  the  variability 
in  the  published  data  on  fatigue  crack  growth  kinetics  might  be 
attributed  to  sources  described  herein.  Caution  is,  therefore, 
recommended  in  attempting  to  draw  statistical  inferences  regard- 
ing material  variability  from  these  data. 
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Table  2:  The  Influence  of  Data  Processing  Procedure 
on  Fatigue  Crack  Growth  Data 


Method  C:  5-point  incremental  polynomial  method. 

Hethoil  D:  Simpson's  rule  (3-point  parabolic  fit). 

Method  E:  Secant  method. 

1 » Same  number  of  data  points  as  Method  A;  2 = all  data  points. 


Table  3;  Statistical  Comparison  of  C and  n Values  from 
Different  Data  Processing  Methods  Using  the 
Same  Number  of  Processed  Data  Pairs  (see  Table  2) 
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I = equal  to,  L = less  than  and  6 = greater  than  on  a 95  percent  confident 
'or  example,  C for  the  5-point  incremental  polynomial  method  (Cl)  is  grea 
;hat  for  the  9-point  incremental  polynomial  method  (Al)  at  the  95  percent 
lence  level.  I = identity. 
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Table  4:  Statistical  Comparison  of  C and  n Values 
from  Different  Data  Processing  Methods 
with  the  Initially  Assumed  Values  (see 
Table  2) . 


Measurement 

Interval 

and 

Precision 

Aa=0. 025  in. , 

a»0.005  in. 

Aa=0 . 025  in. , 

a=0. 00125  in. 

Methods 

C»4xl0"® 

n«2.25 

C»4xl0“^ 

n=2.25 

A1 

L 

G 

L 

G 

B1 

L 

G 

L 

G 

Cl 

L 

G 

L 

G 

D1 

E 

E 

L 

G 

El 

E 

E 

E 

E 

A2 

L 

G 

L 

G 

B2 

L 

G 

L 

G 

C2 

L 

G 

L 

G 

D2 

E 

E 

E 

E 

E2 

E 

E 

E 

E 

See  notes  for  Tables  1,  2 and  3.  L,  E and  G indicate  that  the 
processed  C and  n values  are  less  than,  equal  to  or  greater  than 
the  initially  assumed  values  of  C and  n,  respectively,  at  the  95 
percent  confidence  level. 


( I ) Primary  Data 


(2)  Processed  Data 


■^  = C(AK)" 

dN 


log  AK 

(3)  Empirical  Modeling 


(4)  Prediction 


Figure  1:  Schematic  illustration  of  the  sequential  steps 
involved  in  fatigue  crack  growth  analysis  for 
design. 
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Figure  2:  Illustration  of  procedure  for  assigning  measure' 
ment  error. 
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